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We study the propagation of anisotropic sound and shock waves in dipolar Bose-Einstein condensate in three dimensions 
(3D) as weh as in quasi-two (2D, disk shape) and quasi-one (ID, cigar shape) dimensions using the mean-field approach. 
In 3D, the propagation of sound and shock waves are distinct in directions parallel and perpendicular to dipole axis with 
the appearance of instability above a critical value corresponding to attraction. Similar instability appears in ID and 
not in 2D. The numerical anisotropic Mach angle agrees with theoretical prediction. The numerical sound velocity in 
all cases agrees with that calculated from Bogoliubov theory. A movie of the anisotropic wave propagation in a dipolar 
condensate is made available as supplementary material. 
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1. Introduction 

The alkali metal atoms used in early Bose-Einstein con- 
densate (BEC) experiments have negligible dipole mo- 
ment. However, many bosonic atoms and molecules have 
large dipole moments and ^^Cr [1 and ^^^Dy [2 BECs, 
with a larger long-range dipolar interaction superposed 
on the short-range atomic interaction, have been realized. 
Other atoms like ^^^Er [3 or molecule like ^Li-^^^Cs [4 
with even larger dipole moment are candidates for BEC. 
The superposition of short-range atomic and long-range 
non-local anisotropic dipolar interaction makes the study 
of dipolar BEC (DBEC) very challenging because of the 
appearance of many peculiar properties [1, 5 . 

A tiny object can be dragged in a BEC superfluid below 
the Landau critical velocity [6 without causing any change 
in the superfluid. Above this critical velocity, collective 
excitations are generated in the BEC in the form of soliton 
[7] , turbulent vortex i8j vortex-anti- vortex pair [9 , and 
shock wave [10 , etc. The critical velocity for generation 
of these diverse excitations have different values. It is 0.43 
times the sound velocity for generation of vortices, sound 
velocity for the sound waves, 1.44 times the sound velocity 
for solitons [11 . The original Landau criterion [6 refers 
to the linear excitations only (phonons and rotons in He 

II). 

One remarkable property of a DBEC is its anisotropic 
superfluidity [9 . In a dilute DBEC, or in a dense nondipo- 
lar BEC, like liquid helium in bulk, sound velocity could be 
greater than Landau critical velocity due to roton-like col- 
lective excitations [T?. The Landau critical velocity could 
be anisotropic in a DBEC due to anisotropic rotons [9^. 
Also, anisotropic soliton can be generated in a DBEC [13]. 
Anisotropic collapse has been observed in a ^^Cr DBEC 



[14]. The stability of a trapped DBEC [15 also shows dis- 
tinct anisotropy with a disk-shaped trap leading to more 
stability than a cigar-shaped one [16 . 

Here we consider another manifestation of anisotropic 
superfluidity in a DBEC, e.g., anisotropic sound and shock 
waves. Shock waves have been widely investigated in dif- 
ferent systems, such as, on water surface, in supersonic jet 
and bullet [17J flights, in a gas bubble driven acoustically 
[18 , in a photonic crystal [19 , in the nonlinear Schrodinger 
equation [20] and in BEC of atoms [10], and recently ob- 
served in BEC of polaritons [7]. Using the mean-field 
Gross-Pitaevskii (CP) equation, we study the propagation 
of sound and shock waves in a uniform dilute DBEC in 
three dimensions (3D). Sound propagates anisotropically 
in a 3D DBEC of small dipole moment with larger velocity 
along the axial polarization direction (z) and smaller veloc- 
ity in radial direction (p = {x^y}). Consequently, sound 
wave emitted from a point has a non-spherical ellipsoid- 
like front. We study the anisotropic (oblique) waves when 
a tiny object is dragged along the z and x axes with shock 
and hypersonic velocities. The anisotropic Mach angle [17] 
for drag along z and x axes are in agreement with a theo- 
retical prediction. For a critical dipole moment, the sound 
velocity in x-y plane falls to zero and for larger dipole 
moments an instability due to attraction begins. 

Next we study the sound propagation in trapped DBEC, 
which is of interest in experiments. The effect of dipole mo- 
ment is more prominent in the cigar (ID) and disk (2D) 
shapes with strong traps in radial and axial directions. In 
a spherically- symmetric 3D trap, the effect of dipole mo- 
ment is less pronounced In the cigar shape, there is 
extra attraction due to dipole moment, and the disk shape 
has added repulsion. We study these ID [ST and 2D [22] 
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DBEC using reduced GP equations, where the radial and 
axial variables are integrated out. In the disk-shaped 2D 
DBEC in the radial plane, dipole moment contribute to 
extra repulsion and this makes the system more stable. 
Consequently, sound propagates isotropically for all val- 
ues of dipole moment with a velocity larger than that of 
a non-dipolar EEC. In the cigar-shaped ID DBEC along 
the axial direction, dipole moment contribute to extra at- 
traction and sound propagates with a velocity smaller than 
that of a non-dipolar EEC for dipole moment below a criti- 
cal value. Above this critical value, instability appears due 
to attraction. The numerical sound velocities in ID, 2D 
and 3D DBEC are in agreement with Bogoliubov theory. 

2. Sound and shock waves in Dipolar BEC 

We study a uniform DBEC using the GP equation [23^ 



dt 



2 



dk 



7d(k) 



</>d(r,t), (1) 



where /d(k) = hd(k)Ud(k)^ d = 1,2,3 represent ID, 2D, 
and 3D, respectively, with Ud{k) the short-range plus 
dipolar interaction in momentum space and nrf(k) = 
/ (ire~**^-^|^(^(r)|^ the momentum-space density. The 
space (r) and momentum (k) vectors have d components. 
In 3D, the interaction potential is [23] 

UsCk) = Ana + Airaddji^ ^os^ 



1), 



(2) 



with cos 6 = kz/k, where 6 is the polar angle. Here 
0(^(r,t) is the wave function at time a the atomic scat- 
tering length (taken to be positive here), dipole interaction 
strength add = /ioM^^/(127rfi^), ja the (magnetic) dipole 
moment of an atom, and /io the permeability of free space, 
1>7>— l/2isa tuning parameter controlled by rotating 
orienting fields [22l [24] . In the following we use 7 = 1 and 
shall make some comments on the consequence of negative 
7 on sound propagation. In Eq. ([T]) length is measured in 
units of /o = 1 and time t in units of to = mlo'^/h, 
where m is the mass of an atom. 

The quasi-2D shape is achieved with a strong harmonic 
trap in the axial z direction with oscillator length dz = 
^Jh|muJz^ where ujz is the angular frequency of axial trap. 
In this case, assuming that the axial excitations are in the 
oscillator ground state, the z dependence can be integrated 
out and in Eq. ([T]), we have [22 



/72(k) = 



47ra ^-i^add , ( kpdz 
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dz dz 
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with h2D{0 = 2 — 3v^^e^^erfc(^). The quasi-lD shape is 
achieved with a strong harmonic trap in radial direction 
with oscillator length dp = y/h/mujp, where Up is the an- 
gular frequency of radial trap. Then, the p dependence 
can be integrated out to obtain Eq. ([T]) with |2T] 
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with sioiO = Jo" e-^[3e/{u + e) - Mdu. 



2.1. Bogoliubov spectrum 

The Bogoliubov spectrum of a uniform gas of density Ud 

is mm 



c^d(k) = Viky2) [P/2 + 2/7,(k)n^]. (5) 
In 3D, the Bogoliubov velocity of sound is [9l [T2] 
^3(k) 



c{0) = lim 



IkKo |k| 



coVl + add(3cos2 6> - l)/a, (6) 



where Cq = ^/47^ans is the sound velocity in nondipolar 
medium {add = 0) where ns is density in 3D. This velocity 
is real for 3 cos 19^ > 1. For 3cos(92 < 1 ((9 > 54.7°), the 
3D velocity (|6| can be imaginary for add above a critical 
value, signaling an instability. 

In 2D, the interaction potential (|3| is always positive 
(repulsive) allowing sound propagation in radial direction 
for all add and the Bogoliubov sound velocity is 



lim 

IkKo 



^2(k) 



2n2V27r{a^2add)/dz, (7) 



where 712 is 2D density. If 7 = —1/2 in Eq. ([2|, Cp can be 
imaginary corresponding to the attraction instability. 

In ID, the interaction potential Q, in the kz ^ limit, 
is negative (attractive) for add > ci- The Bogoliubov sound 
velocity is 



lim 

IkKo 



c^i(k) 



^2ni{a - add)/d] 



(8) 



where rii is ID density. The ID sound velocity is imagi- 
nary for add > CL when attraction instability appears. This 
instability resulting from attraction can be avoided by tak- 
ing 7 = -1/2. 

We also calculated Gaussian variational Lagrangians 
[2F for the ID and 2D GP equations [23 and found ef- 
fective interactions Ud{k.), which yielded the variational 
sound velocities identical to the Bogoliubov results ([8| and 

0. 

2.2. Numerical simulations 

We solve the ID, 2D, and 3D GP equations using 
real-time propagation with Crank-Nicolson method |26] 
to study sound and shock waves In 3D, we con- 

sider a generic dipolar atom for numerical simulation with 
a = 4 nm and background density ns = 100 /im~^ = 10^^ 
cm~^. To study sound propagation, a 3D Gaussian pulse 
is placed at the center of the uniform 3D background 
density given by (/)i(r,0) = [100 + 406-^'/^^^')] /im-^, 
w = 2 /im subject to the weak expulsive Gaussian po- 
tential V(r) = 0.00001e-^'/(2^') /im-^ at t = 0. On real- 
time evolution of the GP equation, the 3D Gaussian pulse 
expands. At large times a spherical (nondipolar BEC) 
or an ellipsoid-like (DBEC) sound wave front propagates 
outwards with a central uniform density. From the time 
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Figure 1: Contour plot of sound wave pulse (a) |(/>3(x, 0, z)\'^ {add = 
0), (b) |(/)3(a:,y,0)|2,(add = 30ao), (c) IM^^O, z)\^ , (add = 30ao), 
for sound propagation from center outwards at different times t = 
2, 4, 6, 8, 10. (d) Numerical (num) and Bogoliubov (theory) sound 
velocity for different add- (e) The sound velocity c(6) for different 
polar angle 0. 



evolution of the sound wave front the anisotropic sound 
velocity in different directions is calculated. 

Now we study sound propagation in 3D for a uniform 
DBEC. In Fig. [l](a), we show isotropic sound propagation 
in x-z plane for a nondipolar BEC. In Fig. [l](b), we show 
isotropic sound propagation in x-y plane [9 for a DBEC 
of add = 30ao. In Fig. [l] (c), we show anisotropic sound 
propagation in x-z plane for add = 30ao. The propagation 
along z direction in a DBEC has a velocity larger than that 
for a nondipolar BEC, whereas that in the x-y plane for 
a DBEC has a velocity smaller than that for a nondipolar 
BEC. The Bogoliubov sound velocity for nondipolar atoms 
is 2.242 jam/to compared to the numerical velocity of 2.37 
jiia/tQ. For add = 30ao, the axial Bogoliubov sound ve- 
locity is c(0) = 3.002 fim/to (numerical 3.20 /im/to), and 
the radial Bogoliubov sound velocity is c(7r/2) = 1.741 
/im/to (numerical 1.83 /im/to). In Fig. [l](d), we show the 
axial and radial sound velocities, c(0) and c(7r/2), versus 
add from numerical simulation and Bogoliubov theory with 
good agreement between the two. The radial velocity goes 
to zero for add = and for add > <^ attraction instability 
begins. In Fig. [l](e), we show sound velocity c{0) versus 
for add = 30ao. 
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Potentials 



-2 



/im 



or 

Q 2e-[^'+(^-^.t)2]/o.08 ^^-2 velocity ' and v, 
were used to study shock waves [10 in uniform 3D 
condensates. These potentials can be created by moving 
blue-detuned lasers [9 . These potentials generate waves 



0.0 



0.1 



0.3 

1/M 



0.5 



0.6 



Figure 2: Contour plot of \(j)3(x,0, z)\'^ for wave propagation in a 
uniform 3D BEC for laser drag (a) 6 /im/to and (d) 30 iim/to. The 
same in a DBEC with add = 30ao for drag (b) 6 /im/to and (e) 30 
/im/to along z axis and for drag (c) 6 fim/to and (f) 30 fim/to along 
X axis. The angles a and /3 are shown in (a) and (c). (g) The sina 
versus 1/M plot from theory (|9| and numerical simulation (num) for 
different drag Vx and Vz along x and z axes in the DBEC. 
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Figure 3: Numerical (num) velocity of sound propagation Cz and Cp 
in ID and 2D, respectively, compared with Bogoliubov theory and 
variational (var) approximation. 



in x-z plane along x or z direction. Typical contour plot 
in x-z plane of the 3D isotropic oblique wave is shown 
in Figs. [2] (a) and (d) for a BEC with Vz = Q /im/^o 
(supersonic) and 30 /im/to (hypersonic). Anisotropic 
waves for a DBEC with add = 30ao for Vx^Vz = 6 /im/^o 
and 30 /im/to are displayed in Figs. [2](b), (c), (d), and 
(e). The Mach angle a [17 is related to laser velocity 
v{= Vx^Vz) and polar angle /3 of sound propagation by 



3 



[see, Figs. [2] (a) and (c)] 

sin a = 1/M, M = v/c(/3) . (9) 

In Figs. [2] (a) and (d), M = 6/2.37 = 2.53, and M = 
30/2.37 = 12.66, respectively. For a DBEC, sound veloc- 
ity is different along x and z axes and a will be different 
for laser drag along these axes. For drag along x axis 
c(/3) = c(a), and for drag along z axis c{p) = c(7r/2 — a). 
From a full 3D simulation, for drag along x axis, in Fig. 
[2] (c), a ^ 31° and M v^/c{a) 6/2.99 ^ 2.01 and 
for drag along z axis, in Fig. [2] (b), a ~ 18.3° and 
M = vJc{7r/2 - a) = 6/2.07 ^ 2.90. At supersonic 
velocity, the wave with ripples outside the Mach angle 
is oblique, viz. Figs. [2] (a) — (c). At hypersonic ve- 
locity, the wave, mostly confined in the Mach angle, is 
nearly normal, viz. Figs. |2](d) - (f). To test Eq. ^ 
we did simulation for different Vx{vz) and calculated a 
from 2D contour plots of wave propagation along x and 
z axes for a DBEC with add = 30ao and obtained c(/3) 
and M using results in Fig. [l] (e). In Fig. |2] (g) we 
plot the theoretical and numerical results for sin a ver- 
sus 1/M. A movie of the wave of Figs. |2] (b) and (c) 
is available as supplementary material ^ (also available 
at |http: / /www.youtube.com/watch?v= NR9dB7uZr9M ) . 
Two independent simulations reported in Figs. [T] and [2] 
confirm the anisotropic nature of sound propagation. 

Next we consider sound propagation in trapped medium. 
The effect of dipolar moment is more pronounced in cigar 
and disk shapes and we consider these two cases here. In 
the cigar shape, dipolar moment leads to attraction which 
can eventually lead to collapse for sufficiently large dipole 
moment. In the disk shape the system is always repulsive 
and propagation of sound is allowed for all values of dipole 
moment. The sound wave in cigar and disk shapes are 
studied by the reduced ID and 2D GP equations with uni- 
form density by putting at the center an initial Gaussian 
pulse on top of an expulsive Gaussian potential at t = 
as in 3D. The numerical simulation was done with inter- 
actions ^ and Q using a = 4 nm and ns = 10^^ cm~^ 
and dz = I /im 77,2 = nsV^ and dp = 1 /im, rii = 27in^ 
in 2D and ID, respectively. This will lead to the same 
nondipolar sound velocity in ID, 2D, and 3D. With time 
evolution, the sound wave propagates outwards from the 
center. In Fig. |3] sound velocities calculated for ID and 
2D propagations are plotted versus add together with the 
variational and Bogoliubov sound velocities. Instability 
due to attraction in ID appears for add > from Bogoli- 
ubov theory and variational approximation. In numerical 
simulation the instability appears at a slightly higher value 
of add ^ 80ao ^ 4.23 nm in ID. 

In this study we took the polarization along axial z 
axis. The conclusion about sound velocity in ID and 2D 
will change if we consider the polarization along x axis. 
The h2D of Eq. ^ wih be [9 h2D{kpdz/V2) = -1+ 
3^/7^k'^dze^ erfc(^)/(A/2/Cp) yielding the 2D velocity Cp = 

^J2n2^/2^l{a — add)/dz- This velocity will decrease with 



increasing add and for add > it becomes imaginary cor- 
responding to attraction. This instability due to attrac- 
tion could be avoided by considering 7 = —1/2 in Eq. 
([2|. Also, for polarization along x axis, the ID velocity 
Cz will increase with increasing add avoiding the attraction 
instability. 

3. Summary and Conclusion 

We studied anisotropic sound and shock wave propa- 
gation in DBEC using the mean-field GP equation. In 
a uniform 3D DBEC, sound propagates faster along po- 
larization {z) direction than in radial x-y plane. For 
cidd > there is attraction instability in the radial x-y 
plane. Oblique waves with distinct Mach angles are formed 
when a tiny object is dragged with supersonic velocity in 
a DBEC along or perpendicular to the polarization direc- 
tion. In case of a DBEC, the oblique wave is anisotropic 
and the modified Mach angle relation ^ is found to hold. 
In a quasi-2D disk-shaped DBEC, sound velocity is found 
to increase with add^ whereas in a quasi- ID cigar-shaped 
DBEC sound velocity decreases with add and becomes zero 
for add = and there is attraction instability for add > cl- 
The conditions of attraction instability can be reverted by 
tuning 7 = —1/2 in Eq. ([2|. The attraction instability 
can be implemented by tuning the scattering length by 
a Feshbach resonance [16] and/or dipole moment by ro- 
tating orienting fields [24^ and verified experimentally. In 
all cases the theoretical Bogoliubov sound velocities are in 
agreement with numerical simulation. Experimental con- 
firmation is welcome for the fascinating anisotropic sound 
and shock wave propagation reported here. 
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